Stimulated Raman transitions via multiple atomic levels 
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We consider the stimulated Raman transition between two long-lived states via multiple interme- 
diate states, such as between hyperfine ground states in the alkali-metal atoms. We present a concise 
treatment of the general, multilevel, off-resonant case, and we show how the lightshift emerges nat- 
urally in this approach. We illustrate our results by application to alkali-metal atoms and we make 
specific reference to cesium. We comment on some artifacts, due solely to the geometrical overlap 
of states, which are relevant to existing experiments. 
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I. INTRODUCTION 

The stimulated Raman transition is an extremely pow- 
erful tool for laser manipulation of cold atoms and ions. 
By coupling long-lived states via, but never populat- 
ing, radiative states, experimenters can emulate near- 
ideal two-level quantum systems with no significant de- 
cay [1-3]. This technique has been used to measure sub- 
lincwidth features [4, 5] and to construct atomic interfer- 
ometers which, by exploiting photon recoil, create spa- 
tially separated atomic wave packets which are sensitive 
to gravity [6, 7] or fundamental constants [8, 9]. The 
effective two-level system, which emerges from the Ra- 
man problem, can exhibit behavior such as Rabi flop- 
ping [10, 11], can be used for experiments such as Ram- 
sey interferometry [12, 13], and can provide the qubits 
for quantum information processing [14-16]. Sequences 
of Raman pulses can be used to craft arbitrary superpo- 
sitions in systems with numerous metastable states [17] 
and to prepare such systems in particular states prior to 
coherent manipulation [18]. Raman processes have also 
been used to cool atomic samples to far below the photon 
recoil limit [19-22]. 

Throughout the literature, when the Raman transition 
is discussed, the level structure of the atom is often ap- 
proximated to three levels — two metastable states and 
one intermediate (radiative) state. The Raman problem 
is solved for this prototypical case and then extended, 
without proof, to include the multilevel structure of the 
atom by summing over the various possible routes (see, 
e.g., Ref. [7, §2.1]). Here, by including multiple routes 
from the outset, we confirm that this simple approach is 
correct, show how an expression for the lightshift emerges 
naturally from this treatment, and show that the system 
behaves as a two-level system with an effective coupling 
strength and an effective detuning. 

There is much existing work related to this problem. 
The three-level (single intermediate state) off-resonant 
case has been treated [23] , there have been extensions to 
four levels [24], and the general multilevel problem has 



been recast into "serial" and "parallel" cases [25]. It has 
never been shown rigorously, however, that the three- 
level case can be extended in the way so often assumed. 
In the following, we use the semiclassical approach but 
alternatively one might consider the Jaynes-Cummings 
model [26, 27]. 

This article is structured as follows. We first describe, 
in Sec. II, the Raman transition in a three-level system 
and we show how this can be generalized to include multi- 
ple intermediate states; details of the lengthly calculation 
are confined to the appendix. We then derive, in Sec. Ill, 
expressions for the behavior of the quantum-mechanical 
amplitudes in the general, off-resonance case. In Sec. IV 
we show how these results can be applied to alkali-metal 
atoms, and we conclude in Sec. V. 



II. THREE-LEVEL SYSTEMS 

The simplest system in which a Raman transition may 
be driven is the three-level "A" system, illustrated in 
Fig. 1, in which two long-lived ground states are coupled 
via a radiative upper state which, because the single- 
photon detuning is sufficiently large, is never significantly 
populated. We label the states of the system by |n), with 
states |0) and |2) coupled by the "pump" field of strength 
f2p and frequency wp, and states |1) and |2) coupled by 
the "Stokes" field of strength Q$ and frequency cjg . Using 
the usual correspondence between bra-ket and column 
vector notation (see, e.g., Ref. [28, §II-C]), the Hamilto- 
nian for this system may be represented by the following 
matrix (see, e.g., Ref. [11, §3.2]): 



ujq Qpcosujpt 

U>1 f^scoswst 

fipCOSWpf f^sCOSWs^ 0J2 



(1) 
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Here, as in standard treatments, we assume there is no 
coupling between states |0) and |2) by the Stokes field, 
or between states |1) and |2) by the pump field. 

This Hamiltonian can be simplified by making the ro- 
tating wave approximation and transforming to the inter- 
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FIG. 1: A simple three-level "A" system, in which a Raman 
transition between states |0) and |1) is driven by "pump" and 
"Stokes" fields via intermediate state |2). For each field, the 
coupling strength and frequency are shown in parentheses, 
and frequencies are chosen to be near two-photon resonance: 
up = (uj-2 — ujo) + A; ljs — (cji — ujo) + (A + 8). The single- 
photon detuning A is large compared with the couplings, 
A| f2p,s, and, in this illustration, is negative: A < 0. 
The two-photon detuning 8 is small compared with the sep- 
aration between the ground states and, in this illustration, is 
also negative. 



tion with this Hamiltonian by using unitary transforma- 
tions to find a basis where the time evolution of the states 
is simple and the transformed Hamiltonian is diagonal. 
When one makes such a conversion between bases, is it 
possible to find an equivalent Schrodinger equation with a 
transformed Hamiltonian [11]: if |-0b) = Oba |Va), then 
i(d/dt) \4> B ) = H B \tp B ) where 



O 



BA 



H A 



BA 



l dt BA 



1-1 



(4) 



For the multistatc Hamiltonian Ha in Eq. (3), we 
choose the operator Oba to be the matrix of eigenvec- 
tors; the first term OHO^ 1 is thus the diagonal matrix 
of the eigenvalues. In the appendix, we detail a pro- 
cedure to find the eigensystem of this Hamiltonian; we 
find two eigenvectors which are superpositions of these 
ground states, and N — 2 more which are superpositions 
of the remaining intermediate levels. These N — 2 upper 
states are decoupled from the ground states and so we 
ignore them in the following treatment. 

The difference between the eigenvalues for the two 
ground state eigenvectors is 



action picture. This yields a slowly varying Hamiltonian, 



H 



~i8t 




(2) 



where the pump frequency wp = (uj2 —coq) + A. is detuned 
from single-photon resonance by A, and the difference 
between the pump frequency and the Stokes frequency 
ujs = — + (A + 8) is offset by S from the two- 
photon resonance {uj\ — luq). We now extend this inter- 
action picture to include multiple levels and define the 
Hamiltonian Ha to be 
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The second part of the subscript, 2, 3, . . . , N, denotes the 
level to which the pump or Stokes field couples. Note that 
the oscillation frequency is the same for each Stokes term 
because this depends on the difference in the frequency of 
the fields and not on the Bohr energy of the intermediate 
level. However, the single-photon dctunings do depend 
on the intermediate level Bohr frequencies, but we now 
make the approximation that the detuning is large com- 
pared with the separation of these intermediate levels, 
and hence A := A2 ~ A3 . . . w A^r. 

In this limit, the above Hamiltonian describes a Raman 
system and we expect to see oscillations of population be- 
tween the ground states. We solve the Schrodinger equa- 



n B = 



1 

2A 



P .0*P + i(||n s p-||0 P p) 2 , (5) 



where, for conciseness, we have represented the cou- 
plings as vectors f2p and Sis with components f2p : i and 
rig.j respectively, and have used vector notation for dot 
products and norms. This oscillation frequency is com- 
posed of a coupling strength S1b and a detuning Ab , via 
Hb = \/S! B + A3, analogously to a two-level system: 



B = 
A B = 



|i~2p • rig 1 

2Z\ ' 
l|fts|| 2 -||ftp|| 2 
4A 



(0) 
(7) 



The detuning Ab is readily identified as the lightshift 
and we justify this in the next section. 

The operator which describes the transformation from 
the bare ground states to these dressed ground states can, 
because of normalization, be written as a rotation: 



Oba = 



COS t 

-e~ iSt sin* 



e +l5t sin ( 



(8) 



and the angle 8 is defined by tan (9 = yA B — Ob J /Ob- 

This treatment is sufficient for the on two-photon res- 
onance case, where 5 — 0, but in general the effective 
Hamiltonian H B also contains a time-derivative second 
term, originating from the time-dependence of the oper- 
ator Oba- Away from the two-photon resonance, where 
8 =/= 0, we find the following slowly varying, but never- 
theless time-dependent, effective Hamiltonian: 



H B = 



-6 sin 2 9 -8e +lSt cos 0sin ( 
-Se- lSt cos9sm6 8sin 2 9 + n B 



(9) 
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III. DETUNING FROM RESONANCE 



A. Explicit forms of the amplitudes 



The Hamiltonian Hb in Eq. (9) has the same form as 
that for the simple two-level problem in the interaction 
picture and with the rotating-wavc approximation. We 
can, therefore, use familiar tools to solve this problem. 
First, we transform to find a time- independent Hamilto- 
nian using the operator Ocb in Eq. (4): 



O. 



CB 



1 

e iSt 

He = 



—S sin 2 9 —ScosOsmt 
-S cos 9 sin 9 Scos 2 9 + ili 



(10) 



Next, analogously to the dressed-states approach, we ro- 
tate by an angle 9i (thus defining Odc) where 



The time dependence of the bare state coefficients 
(t) is readily calculable from the procedure described 
above and is stated here for completeness: 



Mt) = A,(o) 



(|n D t) 



A 1 {t)= L4i(0) 
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; — — sin I TjliDT 



fir 



A 1 (0)i£^sin(ifi D i) 
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+ i— — sm I ^iiu'i 
^D 



-iSt 



(15a) 



(15b) 



tan(20 2 ) = 



<5sin(20) 
n B - 5cos{29) 



(11) 



to find a diagonal Hamiltonian Hd. The difference be- 
tween the diagonal elements of Hp, corresponds to the 
phase evolution frequency Qp, of the states in this ba- 
sis. As in the previous section, we see that this oscil- 
lation frequency is composed of a coupling strength Q,p 
and a modified effective detuning Ad: ^d = \A^b + A d , 
where Ad = Ab — 6 is the detuning relative to Ab , which 
was previously identified as the lightshift. 

We now relate the pure phase evolution in this doubly 
dressed basis to the evolution of the bare states by con- 
catenating the transformations that led us to this final 
Hamiltonian: 



O 



DA 



Odc ■ O, 



CB 



Oba 



(12) 



and |^ D ) = D a |Va_), 



cos{ 
sin(' 
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e iSt cos( 
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(13) 

where Dq^i and ^4o,i are the ground- and excited-state 
components of the doubly dressed wave function |V , d) 
and the bare (interaction picture) wave function (V'a)) 
respectively. Finally, the time evolution of the doubly 
dressed states is simply 



D (t) 
Di(t) 



D (t = 0) 
D^t = 0)e lf2r 



(14) 



Using Eq. (13) we can find the dressed state initial 
conditions £*o,i(^ = 0) in terms of the bare state initial 
conditions -Ao,i(f = 0). Using these values, we can then 
use Eq. (14) to find the dressed state coefficients at some 
later time. Finally, we can invert the transformation in 
Eq. (13) to find the time evolution of the bare state am- 
plitudes. 



Hence the system behaves as a two-level system with the 
coupling strength ftp. = and detuning Ad = Ab — S, 
relative to the effective detuning Ab. This justifies our 
previous identification of Ab with the lightshift. 



B. Oscillation amplitude 

The complete, but cumbersome, formulas in Eq. (15) 
describe the behavior of the bare-state amplitudes in 
terms of the bare state initial conditions. If, instead, 
we express this evolution in terms of the initial values 
in the doubly dressed basis, we see clearly that the evo- 
lution is composed of a time-independent offset and an 
oscillation: 



A (t) = cos(0 + 9 2 )D (0) - sin(0 + 2 )£>i(O)e 



(16) 



The population po(t) = \Ao(t)\ in state |0) therefore 
oscillates with pcak-to-pcak amplitude no greater than 
m = sin [2(9 + 62)] which, expressed in terms of the ef- 
fective coupling strength and detuning, is 



fir 



V^B + A D 



(17) 



This envelope function describes a power-broadened 
Lorentzian, centered on the light shifted frequency dif- 
ference between the ground states. This expression rep- 
resents the maximum possible population transfer, and 
any oscillation will be contained within this envelope. 



C. Comments 

A few specific cases are provided here for illustration. 
First, for f2p = (0, 0, tip) and fig = (0, 0, Q$) we recover 
the well-known results for the three level problem. Next 
we note two interesting cases: for ||Op|| = ||^s!l, the 
lightshift Ab is zero. On the other hand, for |i~2p ■ S~2g | = 
0, the Rabi frequency flp is zero. 
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If, as in this last case, the coupling vectors Jlp and l~ig 
are orthogonal, then the transition is not driven. Exam- 
ples include the trivial case where there is no interme- 
diate state to which both ground states are coupled and 
the case where there are states to which both are cou- 
pled, but where these individual coupling strengths sum 
to zero. However, unless the vectors are orthogonal, it 
is possible to adjust the pump and Stokes field strengths 
to ensure the norms of the vectors are equal, and hence 
that the lightshift is zero. 

Our approach relies on the slow time-dependence of the 
interaction-picture Hamiltonian: we require that the sys- 
tem is near two-photon resonance, as previously stated, 
and that there is no coupling of state |0) by the Stokes 
field or of state |1) by the pump field. If present, these 
cross-coupling terms would cause the off-diagonal terms 
in Ha to oscillate in amplitude as well as phase, and the 
treatment in the appendix would no longer be valid. 

For our treatment to be valid, it must therefore be pos- 
sible to identify clearly which field is resonant with which 
transition (see Refs [1, §13.1] and [11, §3.9] for further dis- 
cussion). First, the detuning must be such that no field 
is close to a single-photon resonance: |A| 3> ||J2p.s|| and 
A ± cjio| 3> ||Op,s||- Additionally, the coupling strength 
must be sufficiently small that each ground state can be 
resolved: u>iq ^ ||I2p,s||- If this last condition is violated, 
the system may still appear Ramanlike and exhibit co- 
herent behavior, but it is not described adequately by 
the treatment in this article. 



IV. ALKALI-METAL ATOMS 

We are able to calculate the Rabi frequency and light- 
shift for two states |0) and |1) coupled via a number of 
upper states. A common embodiment of this situation is 
the coupling of two ground hyperfine states via a man- 
ifold of upper hyperfine states in an alkali-metal atom. 
Indeed, there typically exists many such pairs of states, 
but, as ensured by conservation of angular momentum, 
one state is Raman coupled to at most one other state; 
hence, the total system may be treated as a collection of 
independent pairwise couplings. 

As a typical example, consider the Raman transition 
between the ground hyperfine states, via the radiative up- 
per states, in atomic cesium. The pump and Stokes fields, 
both tuned near to the D 2 transition at 852 nm, couple 
states |6 2 S!/2, F = 3) and 1 6 2 S!/ 2 , F = 4), respectively, 
to the 6 2 P 3 / 2 manifold, and have a frequency difference 
near to the hyperfine splitting of 9.2 GHz [29]. The two 
ground states and the intermediate states are 
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6 2 Si/ 2 ;i ? = 3;m F ) , 

6 2 S 1/2 ;F = 4;ra F + q P - q s ) , and 

6 2 P 3/2 ;F = 2,3,4,5;m F + g P ) , 



(18) 



where gp^g = 0, — 1, +1 are the polarizations of the co- 
propagating pump and Stokes fields and correspond to 



linear and left and right circular polarizations, respec- 
tively; nip labels the Zeeman sublevel, corresponding to 
the projection of the total angular momentum F along 
the quantization axis provided by an external magnetic 
field [30]. Linear polarization, in this context, refers 
specifically to the case of the light electric field paral- 
lel to the quantization axis; if these axes are orthogonal, 
the light field interacts with the atom as though it were 
a superposition of left and right circular polarizations. 

The coupling strengths £lp,s depend not only on the 
light intensities but also on the dipole matrix element for 
the transition. Using the Wigner-Eckart theorem [30-32] 
we can split the overlap integral needed to find this dipole 
matrix element and extract from it a purely geometrical 
term G, leaving a term which embodies the other physical 
details of the transition: 

( 6 2 P 3/2 , F', m' F \p.\ 6 2 S 1/2 , F, m F ) = (J\\fi\\J') 

xG(I,J,F,m F ,J',F',m' F ,q) . (19) 

The "reduced" matrix element, denoted by double bars ||, 
depends on many details of the atom, including nuclear 
mass, and is not easily calculated; it can, however, be 
found experimentally from measurements of the upper- 
state lifetime, as described by Loudon [33, Eq. (2.57)] 
and Dcmtrodcr [34]: 



r 



16tt 3 2 J 



1 \(J\W)f 



3e /iA 3 2 J' + 1 
The second part is the product of geometrical terms 



(20) 



G 



-1 



,2F' + .J+I+m F 



F' 



V'(2F' + 1)(2^+1)(2J+1) 

(21) 



F 

-mp 





J' 




){/< 


F 





where the arraylike symbol in parentheses (...) is the 
Wigner 3-j symbol and the similar term in braces {. . .} 
is the Wigner 6-7 symbol [35, §3.3]; both are closely re- 
lated to the Clebsch-Gordan coefficients. This relation 
is described in detail by Edmonds [36]. The two states 
coupled by the Raman interaction are both in 6 2 Si/ 2 so, 
in the calculations that follow, it is only this geometrical 
term which is relevant. 

We imagine the atom in a region of uniform magnetic 
field and consider an experiment where it is possible to 
adjust the frequency difference in order to sweep across 
transitions between various Zeeman sublevels. The prop- 
erties which affect the dipole moment are the various 
quantum numbers: the nuclear spin I = 7/2; the elec- 
tron angular momentum J = 1/2 or J' = 3/2; the total 
angular momentum F = 3 to F = 4 via F 1 = 2,3,4,5; 
and the aforementioned projection mp of F along the 
quantization axis. 

The coupling strength for a dipole transition between 
states \n) and \m) is proportional to the electric field [1]: 



^p,s = -Ep.s {"-I A \ m ) /ft 



(22) 
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where fi is the dipole operator. As above, we can extract 
a geometrical term and, using the vector notation, 



n PiS = E P>s (J\\fi\\J')G P ,s/fi. 



(23) 



Hence, the relative properties of each of the Zeeman sub- 
levels are determined by the geometrical terms Gp,s- 
If we revisit the equations for the coupling strength 
(Eq. (6)) and the lightshift (Eq. (7)), we see that these 
terms appear as |G P ■ Gg| and HGp^H 2 , respectively. 
Thus: 



= {J l^'I? 2 \EpE* s \ |G P • G s | and 



A R = 



2A h 2 

(jm\j') 2 

4A h 2 



£ s | 2 ||G s || 2 -|£p| 2 ||Gp|| 2 



(24) 
(25) 



It is simple to calculate these factors for a given 
initial state and pair of polarizations to examine how 
the lightshift and the coupling strength depend on 
the strength of the applied fields. We find that 
||G|| 2 , for transitions from | J = 1/2, F = I ± 1/2, mp) to 
| </', F', mp + g), driven by light with polarization g, are, 
for an alkali-metal atom with nuclear spin I, given by 



IGI 



1 ± A(J') 



qmp 
21+1 



(26) 



where A(l/2) = -2 for the D 1 transition and 4(3/2) = 1 
for the Di transition. The pump vector Gp couples from 
the lower hyperfine state and corresponds to the nega- 
tive branch (F = I — 1/2); the stokes vector Gs couples 
from the upper hyperfine state and corresponds to the 
positive branch (F = 1 + 1/2). Using these expressions 
and Eq. (25), one may easily calculate the lightshift for 
a given Raman transition in any alkali-metal atom. 

We now turn to the coupling strength fie which, for 
equal polarizations (gp, gs) = (U 1), is symmetrical about 
mp = 0. We find 



IGp-GrI = 



3(27+1) 



(27) 



for the D\ and D2 transitions, where A is given above. 
If we now break this symmetry by choosing, for example, 
(<7p,<7s) = (0, 1), we find 



|Gp • Gs 



3(27+ i; 



y/T(I + 1/2 -m F ). 



(28) 



where T(n) is the n triangular number (1, 3, 6, 10 . . .). 
The dependence of coupling strength on m? level is illus- 
trated in Fig. 2. 

The values for a common arrangement are shown in 
Table I. While of course linear in any overall scaling of 
the intensity, the lightshift has a different dependence 
on the individual field strengths for each of the Zeeman 
sublcvels. It is offset from zero (for unequal intensities) 
and is linear in mp; the lightshift between the hyperfine 
ground states hence has the same dependence on mp as 




FIG. 2: Illustration of the dependence of coupling strength 
described by Eqs (27) (solid) and (28) (dashed). The values 
have physical meaning at integer and half-integer m-p only 
(dots); continuous lines are shown to guide the eye. A line is 
shown for each of several nuclear spins, beginning at / = 1/2 
and increasing in steps of one half out from the center (solid) 
and from bottom to top (dashed). The strongest coupling is 
between extremal mp states by linear and circular polarisa- 
tions. 
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TABLE I: Scaling of the coupling strength and the lightshift 
for the transition |E = 3, tuf) to \F — 4, tuf + qp — qs) for 
(?p, Is) = (1,1), in cesium, in terms of the geometrical parts 
of the dipole matrix elements, as described by Eqs (26) and 
(27). 



the Zeeman shift, and, in, e.g., Rcf. [37], is sufficient to 
account for the majority of the observed spacing between 
the spectral peaks. 

The coupling strengths are not necessarily symmetrical 
about mp = and, for an experiment in which damping 
is important, this may be manifest as a change in the am- 
plitude of the peaks. However, the common arrangement 
of equal polarizations docs not show this asymmetry, and, 
as noted in Rcf. [37] , asymmetry in the preparation of the 
initial states is also important. 



V. CONCLUSIONS 

We have described the stimulated Raman transition 
and, by including multiple intermediate states from 
the outset, have obtained results which give the cou- 
pling strength for the multistate system and from which 
the lightshift naturally emerges. We have applied this 
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method to the cesium atom, given more general expres- 
sions for the alkali-metal atoms, and noted how the lin- 
ear dependence of the lightshift on the rriF level mimics a 
Zeeman shift. We comment on the possibility of the de- 
pendence of the coupling strength onmp level manifest- 
ing as a variation in peak height in experimental spectra. 

Our results were derived in the limit of far detuning 
and our calculations were simplified greatly by this en- 
forced absence of dccoherence. However, for specific cou- 
pling strengths and detunings, the problem of finding 
eigenvalues and vectors can be treated numerically, and 
many efficient algorithms exist for this task. Hence, a 
similar approach might be used for situations including 
coherent population trapping [38, 39] and electromagnet- 
ically induced transparency [40, 41]. 
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Appendix: Finding the eigensystem of the multilevel 
Raman Hamiltonian 



In this appendix, we find the eigenvalues and eigen- 
vectors for the matrix H representing the Hamiltonian 
Ha of our multilevel system, as described in Sec. II. 
For brevity in the derivation, we make the replacements 
x n = §Op, n /fi , y n = \Q s ,ne lSt /fl , and S = A/^o- 
The frequency is, conceptually, the natural frequency 
scale for the problem. 



We first evaluate the term \A 



no 



1. Determinant 

We calculate the determinant of the N x N matrix 
A = H XI, 

-X x 2 x 3 
-A 2/2 2/3 

x* 2 y* 2 S-X .. (A1) 
x\ y% S-X 



to find the characteristic equation and hence the eigen- 
values A of H. If we define A % i to be the matrix A with 
row i and column j removed, and Aij to be the element 
(i,j) of matrix A, then using expansion by minors, 

\A\ = ^(-1)VL„ |A"°| 

= -AL4 00 | + £(-l)"<|A 



iriO 



(A.2) 



1 00 I 



-A y 2 
y*2 S-X 

vt 







2/3 


S-X 



(A 



oo\ 



■£<- 

n>2 



(A' 



Oil 



ra-1,0 



(A.3) 



As before, we decompose the determinant in terms of the 
elements in the first column. The first minor matrix is 
diagonal (A 00 ) 00 = (S — A) J, and hence the first term is 
-X(S-X) N -' 2 . For subsequent terms (A 00 )"" 1 ' !; n > 2, 
we find — j/*y n (<5 — A)^ 3 , where the problem of calcu- 
lating the determinant of each minor matrix is greatly 
simplified by swapping columns to ensure each is upper 
diagonal with diagonal elements {y n ,S — A, 5 — A, . . .). 
Overall, these terms sum to — ||y|| 2 (<5 — X) N ~ 3 . Hence, 



\A™\ = -X(6-\) N - 2 -\\yn6-\) 



\N-3 



We now consider the terms for n > 2: 



X2 x% X4 

-A 2/2 V3 V4 

A 20 = y* S - X 

y% S-X 



(A.4) 



(A.5) 



Before deconstructing A 20 , we note that the next term, 
A 30 , is of the same form if we interchange the columns 



i 30 



and 


x 3 : 















X3 


%4 






-A 


2/2 


2/3 


2/4 






2/2 


S-X 












Hi 








S-X 













X3 


•<2 


•''4 






-A 


2/3 


2/2 


2/4 






2/2 


5 


- A 









Vl 








8 - X 



n>2 



(A.6) 



with a sign change. The first term of this determinant 
has pref actor zero and so can be ignored. For the second 
term, the minor matrix (A 20 ) 10 is upper diagonal and 
hence has determinant |(A 20 ) 10 | = x 2 {S — X) N ~ 3 and, 
similarly, |(A"°) 10 | = (-l) n x n (8 - X) N ~ 3 for n > 2. 

We now investigate the determinants |(A 20 )" l0 | for 
m > 2. In the expansion of |(A 20 ) 20 |, the only terms 
which have nonzero coefficient are | ((A 20 ) 20 ) 00 ] = y 3 {8 — 
X) N - i and |((A 20 ) 20 ) 10 | = ^-A)^- 4 . We can extend 
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this treatment for n > 2 and m > 2. Finally, we find 
\A\ = \ 2 (5-\f- 2 + \(\\x\\ 2 +\\y\\ 2 )(5-\) N -* 

+ \Y.\ x w- x ^\ 2 ( 5 - x ) N ~ A ■ ( A - 7 ) 



'■./ 



2. Eigenvalues 



The equation |A| = clearly has solution A = 5 with 
multiplicity N — 4. With this factor removed, and us- 
ing ||a:|| 2 ||y|| 2 -|a: • y*\ 2 = |Ei,j \ x iVj ~ x jVif to phrase 
this equation in terms of vectors, we obtain 



a 2 (5-a) 2 + a(H 2 + NI 2 )(£-a) 

+ ||:r|| 2 ||y|| 2 -|z-y > f =0, (A.8) 

which is a fourth-order polynomial in A (with leading 
coefficient 1) and hence the product of the remaining 
roots Xi equals the constant term. This term is finite, 
and so at least one A^ becomes negligible as 6 — > ±oo, 
and hence we can make the approximation (6 — A) — > 6 
in this limit. The resulting equation is a quadratic in A 
with solutions 



-(ll.rf + llyll^i^dl^P-llyll^+^.y*! 2 

A± ~ 26 ■ 

(A.9) 

We now seek the eigenvectors associated with these 
two finite eigenvalues. 



3. Eigenvectors 

The eigenvalue equation Ha = Xa yields the following: 
^^.TiOi = Aao; (A. 10) 



i>2 



'^2yiai = Aai; 



(A.ll) 



i>2 



x*a + y*ai = (A - S)a l for i > 2 . (A. 12) 

If we multiply Eq. (A. 12) by a* and sum over i > 2, and 
then enforce the normalization condition E a*ai = 1, we 
arrive at 

(A -S) [1 - {\a \ 2 + M 2 )] - A (|a | 2 + M 2 ) , (A.13) 
and hence 

|a | 2 + |ai| 2 = l-A/5, (A.14) 

which, in the limit of large 6, tends to unity. By way of 
confirmation, we see from Eq. (A. 12) that 



x*a + y*ai . 
ai = for i > 2 . 

A — 6 



(A.15) 

which tend to zero in this limit. The two eigenstates 
associated with the two finite eigenvalues are hence or- 
thogonal superpositions of the two ground eigenstates; 



we represent this transformation as a rotation and pro- 
ceed to find its angle. Using Eqs (A. 10) and (A. 12), wc 
have 



a 



i>2 



i 



A(A - 6) 

and similarly for a\. Hence we obtain 
ai a Q y-x* + a 1 y-y 



a x ■ x* + aix ■ y*\ , (A. 16) 



(A.17) 



which, because |a | 2 + |ai| 2 = 1, we can express as the 



do aox ■ x* + a\X ■ y* 
1 |ao| 2 

tangent of an angle: 



e v tanfl = ai/a = — , (A. 18) 

2x ■ y* 



where \ = y 4|cc • y*| 2 + (j|y|| 2 — ||a;|| 2 ) and <j) ~ 
-arg(cc • y*). 

The transformation from the bare-state basis to this 
dressed-state basis can hence be described by the rotation 



1+) 
h) 



cos 6» e+^sin(9 \ / |0) \ , 
-e-^srne cos9 ) I |1) ) ' [ } 



and oscillations are thus driven with amplitude m = 
2 cos 9 sin 9 at the rate 

O B = O (A + -A_) = ^| x . (A.20) 

We identify the effective coupling strength Ob and de- 
tuning Ab in terms of the angle 9 defined above: 



Ob = sin 29 Ob and 
Ab = cos 29 B . 



(A.21) 
(A.22) 



l-cos 29 



Using the trigonometric identity tan 9 - 
obtain 



wc 



and, finally, 



sin26» = 2\x-y*\/x and (A.23) 
COS20 = {\\y\\ 2 -\\x\\ 2 )/ X , (A.24) 



Q2 

B = -^2|a:-y* and (A.25) 
Ab = f (j|y|| 2 -M| 2 ). (A.26) 
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